Quasi-Degenerate Neutrinos with non-zero #13 in 
Type I+II Seesaw Models 

Debasish Borah* and Mrinal Kumar Das^ 
Department of Physics, Tezpur University, Tezpur-784028, India 

Abstract 

We study the survivability of quasi-degenerate neutrino (QDN) masses with normal as well as 
inverted hierarchical mass patterns in the presence of both type I and type II seesaw contributions 
to neutrino mass in generic left-right symmetric models. At leading order, the Dirac neutrino mass 
matrix is assumed to be diagonal with either charged lepton (CL) type or up quark (UQ) type 
structure which gets corrected by non-leading effects giving rise to deviations from tri-bi-maximal 
1 (TBM) mixing and hence non-zero value of #13. Using the standard form of neutrino mass matrix 

; which incorporates such non- leading effects within flavor symmetry models like A4, we parametrize 

the neutrino mass matrix incorporating both oscillation as well as cosmology data. Also considering 
extremal values of Major ana CP phases such that the neutrino mass eigenvalues have the structure 
(mi, —777-2,7713) and (7771,7772,7773), we then calculate the predictions for neutrino parameters in the 
presence of both type I and type II seesaw contributions, taking one of them dominant and the 
other sub-dominant. We show that QDN scenario can still survive in our framework with certain 
exceptions. 
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I. INTRODUCTION 



Recent neutrino oscillation experiments have provided significant amount of evidence 
which confirms the existence of the non-zero yet tiny neutrino masses [1] . We know that the 
smallness of three Standard Model neutrino masses can be naturally explained via seesaw 
mechanism. In general, such seesaw mechanism can be of three types : type I [2], type II 
[3] and type III [4]. All these mechanisms involve the inclusion of additional fermionic or 
scalar fields to generate tiny neutrino masses at tree level. Although these seesaw models 
can naturally explain the smallness of neutrino mass compared to the electroweak scale, we 
are still far away from understanding the origin of neutrino mass hierarchies as suggested by 
experiments. Recent neutrino oscillation experiments T2K [5], Double ChooZ [6], Daya-Bay 
[7] and RENO [8] have not only made the earlier predictions for neutrino parameters more 
precise, but also predicted non-zero value of the reactor mixing angle #13. The latest global 
fit value for 3a range of neutrino oscillation parameters [9] are as follows: 

Am 2 21 = (7.00 - 8.09) x 10~ 5 eV 2 

Am 2 x (NH) = (2.27 - 2.69) x 10~ 3 eV 2 
Am 2 3 (IH) = (2.24 - 2.65) x 10~ 3 eV 2 

sin 2 # 12 = 0.27-0.34 

sin 2 # 23 = 0.34 - 0.67 

sin 2 # 13 = 0.016 - 0.030 (1) 

where NH and IH refers to normal and inverted hierarchy respectively. 

The above recent data have positive evidence for non-zero #13 as well, which was earlier 
thought to be zero or negligibly small. Non-zero 6*13 can be explained by incorporating 
various corrections to the standard TBM mixing. The standard TBM mixing pattern is 
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which predicts sin 2 ^ = |, sin 2 ^ = \ and sin 2 #i3 = 0. However, since the latest data have 
ruled out sin 2 6 , 13 = 0, there arises the need to go beyond the TBM framework. Since the 



experimental value of #13 is still much smaller than the other two mixing angles, TBM can 
still be a valid approximation and the non-zero #13 can be accounted for by incorporating 
non-leading contributions to TBM coming from charged lepton mass diagonalization, for 
example. There have already been a great deal of activities in this context [10, 11] and the 
latest data can be successfully predicted within the framework of several interesting models. 
These frameworks which predict non-zero #13 may also shed light on the Dirac CP violating 
phase which is still unknown (and could have remained unknown if # 13 were exactly zero). 

Apart from predicting the correct neutrino oscillation data as well as the Dirac CP phase, 
the nature of neutrino mass hierarchy is also an important yet unresolved issue. Use of 
specific grand unified models explaining the seesaw mechanisms has also been done in the 
last few years to study the hierarchy of neutrino masses. An analysis done in [12] showed 
that every normal neutrino mass hierarchy solution of a grand unified model corresponds 
to an inverted hierarchy solution. It was also mentioned in their work that any future 
observation of inverted hierarchy would tend to disfavor the grand unified models based on 
the conventional type I seesaw mechanism. But models with type II and type III or models 
based on conserved L e — — L u symmetry may favor the inverted hierarchical nature of 
neutrino masses. Models based on seesaw mechanism with three right handed neutrinos can 
also generate inverted hierarchical neutrino masses [13] within the framework of bi-maximal 
mixing. Understanding the correct nature of hierarchy can also have non-trivial relevance 
in leptogenesis as well as cosmology. For example, the latest cosmology data on the sum 
of absolute neutrino masses [14] have already ruled out the scenario of quasi-degenerate 
(QDN) neutrino masses with m, > 0.1 eV. From supernova neutrinos point of view, it 
was shown [15] that one can discriminate the inverted hierarchy from the normal one if 
sin 2 #i3 > a few x 10 -4 . If a particular neutrino mass hierarchy is assumed this can bias 
cosmological parameter constraints [16] like dark energy equation of state parameter as well 
as the sum of the neutrino masses. Therefore, the study of normal and inverted hierarchy 
using different types of seesaw mechanism is very important both from neutrino physics and 
cosmology point of view. 

In view of above, the present work is planned to carry out a study of QDN scenario 
with normal and inverted hierarchical neutrino masses in the framework of generic left-right 
symmetric models (LRSM) [17]. This is an extension of our earlier work [18] where we 
considered TBM type /x — r symmetric neutrino mass matrix. In that study, the dominating 
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type I seesaw term was kept fixed and type II term was varied and the predictions for neutrino 
oscillation parameters were calculated. In this present work, we extend the earlier work to 
include non-zero #13 as well as to consider the case where type II seesaw term dominates over 
type I. We parametrize the neutrino mass matrix for QDN scenario using global fit neutrino 
oscillation as well as cosmology data. The dominating seesaw term is then used to find out 
the right-handed Majorana neutrino mass matrix by assuming two different types of Dirac 
neutrino mass matrices (CL and UQ type) and two different Majorana CP phase patterns 
((mi, -m2, 7713) and (mi, m>2, tus)). Fixing the dominating seesaw term in this way, the other 
seesaw term is allowed to vary and predictions are calculated for all these variations. By 
calculating the predictions for neutrino oscillation parameters, we show that both type I 
dominating and type II dominating cases give almost identical results. We find that normal 
hierarchy along with CL type Dirac neutrino mass matrix gives rise to predictions within 
3cr range of experimental data for both types of Majorana CP phase patterns. Whereas, 
inverted hierarchy with CL type Dirac mass matrix can exist only with the Majorana CP 
phase pattern (mi, —7712,7713). Also, we show that UQ type Dirac neutrino mass matrix is 
disfavored for both the types of hierarchies and Majorana CP phases. 

This paper is organized as follows: in section II we discuss the methodology of type II 
seesaw mechanism in generic LRSM. In section III we discuss our numerical analysis and 
results. We then finally conclude in section IV. 

II. TYPE II SEESAW IN LRSM 

Type I seesaw framework is the simplest mechanism for generating tiny neutrino masses 
and mixing. Such a mechanism is possible in extensions of the standard model by three 
right handed neutrinos. There is also another type of non-canonical seesaw formula (known 
as type-II seesaw formula) [3] where a left-handed Higgs triplet Al picks up a vacuum ex- 
pectation value (vev). This is possible both in the minimal extension of the standard model 
by Al or in other well motivated extensions like left-right symmetric models (LRSM) [17]. 
The seesaw formula in LRSM can be written as 

m LL = m^ L + m{ L (3) 
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where the usual type I seesaw formula is given by the expression, 



m 
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Here ttilr is the Dirac neutrino mass matrix. The above seesaw formula with both type I 
and type II contributions can naturally arise in extension of standard model with three right 
handed neutrinos and one copy of A^. However, we will use this formula in the framework 
of LRSM where M RR arises naturally as a result of parity breaking at high energy and 
both the type I and type II terms can be written in terms of M RR . In LRSM with Higgs 
triplets, M RR can be expressed as M RR = v R f R with v R being the vev of the right handed 
triplet Higgs field A R imparting Majorana masses to the right-handed neutrinos and f R is 
the corresponding Yukawa coupling. The first term mj^ L in equation (3) is due to the vev 
of SU{2) L Higgs triplet. Thus, m 1 ^ = /lVl and M RR = f R v R , where vl jR denote the vev's 
and Jl,r are symmetric 3x3 matrices. The left-right symmetry demands Jr = Jl = f ■ The 
induced vev for the left-handed triplet vl can be shown for generic LRSM to be 



with M\y ~ 80.4 GeV being the weak boson mass such that 

\v L \ « M w « \v R \ 

In general 7 is a function of various couplings in the scalar potential of generic LRSM and 
without any fine tuning 7 is expected to be of the order unity (7 ~ 1). type-II seesaw 
formula in equation (3) can now be expressed as 



With above seesaw formula (5), the neutrino mass matrices are constructed by considering 
contributions from both type I and type II terms. The choice of v R however, remains 
ambiguous in the literature where different choices of v R are made according to convenience 
[19-22]. However, in this present work we will always take v R as v R = 7-^- — 7 x 10 15 GeV 
[22]. It is worth mentioning that, here SU(2) R x U(1)b-l gauge symmetry breaking scale 
(as in generic LRSM) v R is the same as the scale of parity breaking [19]. Using this form of 
v R , the seesaw formula (5) becomes 



vl = 7 




m LL = -i(M w /v R ) 2 M RR - m LR M RR m T LR 



(5) 
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Quantitatively, either of the two terms on the right hand side of equation (6) can be 
dominant or both the terms can be equally dominant. However, for generic choices of 
symmetry breaking scales (mentioned above) as well as the Dirac neutrino mass matrices 
(generically to be of same order as corresponding charged lepton masses), both type I and 
type II term can be equally dominant only when the dimensionless parameter 7 is fine tuned 
to be very small. We check this by equating both the terms to the best fit value of ttill as 
we will show in a tabular form II, III in our section on numerical analysis. We skip such a 
fine-tuned case here and consider two other possible cases in our work: one in which type 
I term dominates whereas type II term is present as a small perturbation and the other in 
which type II term dominates with type I term as a small perturbation. 

A. Case I: Dominating Type I Seesaw 

In this case, the second term in the equation (6) gives the leading contribution to 
and hence we compute the right-handed neutrino mass matrix Mrr by using the inverse 
type I seesaw formula Mrr = m^ R m~ L \mi J R where we use the best fit ttill and generic 
as will be shown in the section III. Here we hold M&r fixed, so the first term in equation 
(6) is dependent on the value of 7 while second term is fixed. For 7 ~ 1, the first term 
has minimum contribution whereas for smaller values of 7, the contribution of the type II 
term will increase. We vary the dimensionless parameter 7 from 0.001 to 1.0 and check 
the survivability of neutrino mass models with contributions from type I and type II terms. 
We adopt a natural selection for the survival of neutrino mass models which have the least 
deviation of 7 from unity. Nearer the value of 7 to one, better the chance for the survival 
of the model in question. Thus the value of 7 is an important parameter for the proposed 
natural selection of the neutrino mass models in question. 

B. Case II: Dominating Type II seesaw 

This is the case where the first term in the equation (6) gives the leading contribution to 
m^i. This scenario within grand unified models like S'O(IO) have been discussed in [23]. In 
this case, we compute M^/7 from the inverse type II seesaw formula 




(7) 
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TABLE I: Input parameters and predictions for parameters in the neutrino mass matrix 



Parameters 


IH(+-+) 


IH(+++) 


NH(+-+) 


NH(+++) 


Am^lO" 5 eV 2 ] 


7.50 


7.50 


7.50 


7.50 


lAm^llKT 3 eV 2 ] 


2.43 


2.43 


2.47 


2.47 


sm 2 6*i2 


0.340 


0.340 


0.340 


0.340 


sin 2 6*23 


0.606 


0.606 


0.606 


0.606 


sin 2 0i3 


0.022 


0.022 


0.022 


0.022 


m 3 (eV) 


0.075 


0.075 


0.090 


0.101 


m2 (eV) 


-0.090 


0.090 


-0.075 


0.089 


mi (eV) 


0.089 


0.089 


0.075 


0.088 


Et m i ( eV ) 


0.25 


0.25 


0.24 


0.28 


X 


0.0206 


0.0803 


0.0251 


0.0893 


y 


-0.0551 


0.0045 


-0.0503 


-5.3969 x 10~ 6 


z 


0.0068 


0.0068 


0.0069 


0.0060 


w 


0.0015 


0.0015 


0.0015 


0.0013 



It should be noted that here we are keeping Mrr/^ constant instead of just Mrr as in 
the case I. Thus, although the first term in equation (6) remains constant, the second term 
varies as we vary 7 and have the minimal contribution for 7 ~ 1. Similar to the case I, here 
also we vary 7 between 0.001 and 1.0 and calculate the predictions for neutrino parameters. 
The details of both these case will be presented in details in the next section. 



III. NUMERICAL ANALYSIS AND RESULTS 



For the purpose of numerical analysis we use the following parametrization of neutrino 
mass matrix which is a combination of TBM version of \i — r symmetric matrix with the 
addition of correction terms coming from charge lepton mass matrix diagonalization [11]. 

x y — w y + w ^ 
f^LL =y — wx + z + wy — z (8) 
y + w y — z x + z — 



\ 



J 
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TABLE II: Right Handed Majorana Neutrino masses in GeV for type I seesaw dominance 



(m, n) 


(6,2) 


(8,4) 








( 7.35 x 10 7 


2.40 x 10 s 


4.79 x 10 10 ^ 










f 1.87 x 10 8 


1.96 x 10 8 


1.22 x 10 11 ^ 






IH(+++) 






2.40 x 10 s 
^ 4.79 x 10 10 


7.10 x 10 10 
9.10 x 10 10 


9.10 x 10 19 
3.14 x 10 13 j 










1.96 x 10 8 
^ 1.22 x 10 11 


6.25 x 10 8 
1.27 x 10 11 


1.27 x 10 11 
8.00 x 10 13 j 










/ 


1.47 x 10 7 


-1.46 x 10 9 


1.36 x 10 10 


\ 






/ 


6.71 x 10 7 


-1.06 x 10 8 


4.43 x 10 10 


\ 




IH(+-+) 






-1.46 x 10 9 


2.13 x 10 10 


-9.04 x 10 11 








-1.06 x 10 s 


-1.35 x 10 8 


-6.87 x 10 10 








V 


1.36 x 10 10 


-9.04 x 10 11 


1.15 x 10 13 


J 






I 


4.43 x 10 10 


-6.87 x 10 10 


2.92 x 10 13 


1 










' 5.27 x 10 7 


2.58 x 10 s 


3.86 x 10 10 N 










' 1.32 x 10 8 


1.47 x 10 8 


9.81 x 10 10 N 






NH(+++) 






2.58 x 10 8 
^ 3.86 x 10 10 


6.54 x 10 10 
1.24 x 10 11 


1.24 x 10 11 
2.86 x 10 13 j 










1.47 x 10 8 
^ 9.81 x 10 10 


5.49 x 10 8 
1.08 x 10 11 


1.08 x 10 11 
7.27 x 10 13 j 










/ 


9.86 x 10 6 


-1.65 x 10 9 


1.07 x 10 10 


\ 








5.82 x 10 7 


-1.35 x 10 8 


3.79 x 10 10 


\ 




NH(+-+) 






-1.65 x 10 9 


1.66 x 10 10 


-9.95 x 10 11 








-1.35 x 10 s 


-2.09 x 10 8 


-8.60 x 10 10 








I 


1.07 x 10 10 


-9.95 x 10 11 


9.69 x 10 12 


) 






I 


3.79 x 10 10 


-8.60 x 10 10 


2.46 x 10 13 


) 





TABLE III: Mrr/j in GeV for Type II seesaw dominance 



IH(+++) 




' 2.27 x 10 16 8.68 x 10 14 1.73 x 10 15 N 

8.68 x 10 14 2.51 x 10 16 -6.40 x 10 14 
^ 1.73 x 10 15 -6.40 x 10 14 2.43 x 10 16 y 




IH(+-+) 




^ 5.86 x 10 15 -1.60 x 10 16 -1.51 x 10 16 N 

-1.60 x 10 16 8.24 x 10 15 -1.75 x 10 25 
i -1.51 x 10 16 -1.75 x 10 16 7.37 xlO 15 




NH(+++) 




' 1.26 x 10 16 -1.92 x 10 14 1.90 x 10 14 
-1.92 x 10 14 1.37 x 10 16 -8.57 x 10 14 
i 1.90 x 10 14 -8.57 x 10 14 1.33 x 10 16 




NH(+-+) 




' 7.14 x 10 15 -1.47 x 10 16 -1.38 x 10 16 
-1.47 x 10 16 9.56 x 10 15 -1.62 x 10 16 
^ -1.38 x 10 16 -1.62 x 10 16 8.67 x 10 15 J 
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TABLE IV: Predictions for neutrino parameters for CL type mtR with IH(+-+) 



Parameters 


Case I 
7 = 0.25 


Case II 
7 = 0.25 


Case I 
7 = 0.50 


Case II 
7 = 0.50 


Case I 
7 = 0.75 


Case II 
7 = 0.75 


Case I 
7 = 1.00 


Case II 
7 = 1.00 


Am^[l(T 5 eV 2 ] 


9.84 


9.72 


8.64 


8.59 


8.26 


8.22 


8.08 


8.05 


Am^lCT 3 eV 2 ] 


2.43 


2.43 


2.43 


2.43 


2.43 


2.43 


2.43 


2.43 


sin 2 6*23 


0.603 


0.604 


0.605 


0.605 


0.605 


0.605 


0.605 


0.606 


sin 2 6»i2 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


sin 2 6»i3 


0.021 


0.021 


0.022 


0.022 


0.022 


0.022 


0.022 


0.022 



TABLE V: Predictions for neutrino parameters for CL type ttilr with IH(+++) 



Parameters 


Case I 
7 = 0.25 


Case II 
7 = 0.25 


Case I 
7 = 0.50 


Case II 
7 = 0.50 


Case I 
7 = 0.75 


Case II 
7 = 0.75 


Case I 
7 = 1.00 


Case II 
7 = 1.00 


Am^JlO" 5 eV 2 ] 


10.06 


9.91 


8.27 


8.23 


7.89 


7.87 


7.75 


7.74 


Am 2 3 [10- 3 eV 2 ] 


2.45 


2.45 


2.43 


2.43 


2.43 


2.43 


2.43 


2.43 


sin 2 (9 23 


0.319 


0.328 


0.446 


0.453 


0.498 


0.502 


0.523 


0.526 


sin 2 6»i2 


0.270 


0.272 


0.306 


0.308 


0.319 


0.320 


0.324 


0.325 


sin 2 6»i3 


0.102 


0.099 


0.062 


0.060 


0.048 


0.047 


0.041 


0.040 



where w denotes the deviation of tull from that within TBM frameworks and setting it to 
zero, the above matrix boils down to the familiar /i — r symmetric matrix. The eigenvalues 
of this matrix are m\ = x + 2y, m 2 = x — y + z — \/ z 2 + 3w 2 and m 3 = x — y + z + V z 2 + 3w 2 
Then we parameterize the above matrix for QDN case. From presently available cosmo- 
logical constraints, the upper bound on sum of neutrino masses has come down to the lowest 
value Xli m « — 0-28 eV [14] which has ruled out QDN neutrino models with > 0.1 eV. 
Parametrization of the matrix (8) is done with this upper bound and taking the largest 
allowed value m,; < 0.1 eV consistent with the latest cosmological data. A classification for 
three-fold QDN neutrino masses [24] with maximum Majorana CP violating phase in their 
eigenvalues is used here. CP phase patterns in the mass eigenvalues for both NH and IH 

are taken as: (m±, — m 2 ,m 3 ) (denoted as (H h)) and (m 1 ,m 2 ,m 3 ) (denoted as (+ + +)). 

Using the best fit values of the global neutrino oscillation observational data [9] on solar 
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TABLE VI: Predictions for neutrino parameters for UQ type mm with IH(+-+) 



Parameters 


Case I 
7 = 0.25 


Case II 
7 = 0.25 


Case I 
7 = 0.50 


Case II 
7 = 0.50 


Case I 
7 = 0.75 


Case II 
7 = 0.75 


Case I 
7 = 1.00 


Case II 
7 = 1.00 


AmlJlO- 5 eV 2 ] 


13.28 


12.99 


10.33 


10.19 


9.37 


9.28 


8.93 


8.86 


Am^lO" 3 eV 2 ] 


2.44 


2.44 


2.43 


2.43 


2.43 


2.43 


2.43 


2.43 


sin 2 6*23 


0.600 


0.600 


0.603 


0.603 


0.604 


0.604 


0.605 


0.605 


sin 2 6»i2 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


sin 2 6»i3 


0.020 


0.020 


0.021 


0.021 


0.021 


0.021 


0.021 


0.021 



TABLE VII: Predictions for neutrino parameters for UQ type with IH(+++) 



Parameters 


Case I 
7 = 0.25 


Case II 
7 = 0.25 


Case I 
7 = 0.50 


Case II 
7 = 0.50 


Case I 
7 = 0.75 


Case II 
7 = 0.75 


Case I 
7 = 1.00 


Case II 
7 = 1.00 


Am^JlO" 5 eV 2 ] 


18.27 


17.70 


11.20 


10.99 


9.39 


9.28 


8.70 


8.63 


Am 2 3 [10- 3 eV 2 ] 


2.51 


2.51 


2.46 


2.46 


2.44 


2.44 


2.44 


2.44 


sin 2 (9 23 


0.148 


0.155 


0.274 


0.283 


0.356 


0.365 


0.406 


0.413 


sin 2 6»i2 


0.200 


0.204 


0.253 


0.256 


0.281 


0.283 


0.295 


0.297 


sin 2 6»i3 


0.157 


0.155 


0.117 


0.114 


0.090 


0.087 


0.074 


0.072 



and atmospheric neutrino mass squared differences, and taking m,; < 0.1 eV, predictions for 
parameters in the matrix (8) are calculated as shown in table I. 

After parameterizing the neutrino mass matrix using oscillation and cosmology data, we 
consider the two cases I and II mentioned in the previous section one by one. First, we 
consider the case I i.e., type I dominance and calculate the right-handed Majorana neutrino 
matrix Mrr using the inverse type I seesaw formula 

M RR = m T LR mllm LR , (9) 

To calculate the Mrr for each case, we need to have the Dirac neutrino mass matrix {ttilr)- 
We take the Dirac neutrino mass matrix ttilr to be diagonal at leading order (LO) with 
either charged lepton mass structure up quark mass structure. The general form of Dirac 
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TABLE VIII: Predictions for neutrino parameters for CL type mm with NH(+-+) 



Parameters 


Case I 
7 = 0.25 


Case II 
7 = 0.25 


Case I 
7 = 0.50 


Case II 
7 = 0.50 


Case I 
7 = 0.75 


Case II 
7 = 0.75 


Case I 
7 = 1.00 


Case II 
7 = 1.00 


AmlJlO- 5 eV 2 ] 


6.35 


6.41 


6.93 


6.96 


7.12 


7.14 


7.21 


7.23 


Am^lO" 3 eV 2 ] 


2.47 


2.47 


2.47 


2.47 


2.47 


2.47 


2.47 


2.47 


sin 2 6*23 


0.604 


0.604 


0.605 


0.605 


0.605 


0.605 


0.606 


0.606 


sin 2 6»i2 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


sin 2 6»i3 


0.021 


0.021 


0.022 


0.022 


0.022 


0.022 


0.022 


0.022 



TABLE IX: Predictions for neutrino parameters for CL type itilr with NH(+++) 



Parameters 


Case I 
7 = 0.25 


Case II 
7 = 0.25 


Case I 
7 = 0.50 


Case II 
7 = 0.50 


Case I 
7 = 0.75 


Case II 
7 = 0.70 


Case I 
7 = 1.00 


Case II 
7 = 1.00 


Am^JlO" 5 eV 2 ] 


9.24 


9.11 


8.06 


8.01 


7.79 


7.77 


7.70 


7.68 


Am 2 3 [10- 3 eV 2 ] 


2.48 


2.48 


2.47 


2.47 


2.47 


2.47 


2.47 


2.47 


sin 2 (9 23 


0.590 


0.590 


0.598 


0.598 


0.601 


0.601 


0.602 


0.602 


sin 2 6»i2 


0.118 


0.124 


0.209 


0.213 


0.249 


0.252 


0.270 


0.271 


sin 2 6»i3 


0.021 


0.022 


0.022 


0.022 


0.022 


0.022 


0.022 


0.022 



HI 



(0) 
LB, 



( 



\ 



) 



m 



f 



(10) 



neutrino mass at LO is 

/ X m 
A n 
1 

where mf corresponds to m T tan / 9 for (m,n) = (6,2), tan/3 = 40 in case of charged lepton 
and m t for (m,n) = (8,4) in the case of up-quarks [25, 26]. A = 0.22 is the standard 
Wolfenstein parameter. Just like the charged lepton mass matrix (which is diagonal at LO) 
gets corrected by non-leading terms in generic flavor symmetric models discussed in [11], 
we also consider similar non-leading contributions to the Dirac neutrino mass matrix. Such 
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TABLE X: Predictions for neutrino parameters for UQ type tulr with NH(+-+) 



Parameters 


Case I 
7 = 0.25 


Case II 
7 = 0.25 


Case I 
7 = 0.50 


Case II 
7 = 0.50 


Case I 
7 = 0.75 


Case II 
7 = 0.70 


Case I 
7 = 1.00 


Case II 
7 = 1.00 


Am^[l(T 5 eV 2 ] 


4.05 


4.24 


5.80 


5.89 


6.37 


6.43 


6.64 


6.68 


Am^lCT 3 eV 2 ] 


2.48 


2.48 


2.47 


2.47 


2.47 


2.47 


2.47 


2.47 


sin 2 6*23 


0.601 


0.601 


0.603 


0.604 


0.604 


0.604 


0.605 


0.605 


sin 2 6»i2 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


0.340 


sin 2 6»i3 


0.020 


0.020 


0.021 


0.021 


0.021 


0.021 


0.022 


0.022 



TABLE XI: Predictions for neutrino parameters for UQ type itilr with NH(+++) 



Parameters 


Case I 
7 = 0.25 


Case II 
7 = 0.25 


Case I 
7 = 0.50 


Case II 
7 = 0.50 


Case I 
7 = 0.75 


Case II 
7 = 0.70 


Case I 
7 = 1.00 


Case II 
7 = 1.00 


Am^JlO" 5 eV 2 ] 


14.53 


14.11 


9.93 


9.76 


8.76 


8.67 


8.32 


8.26 


Am 2 3 [10- 3 eV 2 ] 


2.53 


2.53 


2.49 


2.49 


2.48 


2.48 


2.48 


2.48 


sin 2 (9 23 


0.563 


0.565 


0.585 


0.586 


0.592 


0.593 


0.596 


0.596 


sin 2 6»i2 


0.021 


0.023 


0.090 


0.094 


0.145 
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0.180 


0.185 


sin 2 6»i3 


0.021 


0.021 


0.021 


0.021 
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0.022 


0.022 


0.022 



non-leading contribution can be written as 

/ 



5m 



(i) 

LR 



\ 



(11) 



0(A m ) 0(X m ) G(X m ] 
0(X n ) 0{X n ) 0(X n ) 

\ 0(1) 0(1) 0(1) J 

where £ parametrizes the departure from LO approximation. Since the parameter w in the 
neutrino mass matrix (8) parametrizes the deviation from LO approximation (which is TBM 
for generic A 4 flavor symmetric models), we assume £ ~ w for our calculation. Now, using 
m LR, — m LR + ^ m LR m ^he inverse type I seesaw formula above, we calculate Mrr which 
are presented in the Table II. Using the values of Mjyj in the equation (5), we check the 
variations of mass squared differences as well as the mixing angles with respect to 7 (varying 
from 0.001 to 1). 
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FIG. 1: Variation of the predicted values of Am^ as a function of 7 in NH case 

For case II i.e., type II dominating seesaw, we calculate M RR /^ first using equation (7). 
The values of Mrr/j are shown in table III. This fixes the first term in the seesaw formula 
(6). Since Mrr is not fixed in this case and varry with 7 as 

Axl0 15 \ 2 

the second term in the seesaw formula (6) also varies with 7 with minimum contribution for 
7 ~ 1. Similar to the case I, here also we vary this dimensionless parameter 7 from 0.001 to 
1.0 and check the variations of all the neutrino oscillation parameters. 

The predictions for neutrino oscillation parameters at different values of 7 = 
0.25, 0.50, 0.75, 1.00 for both case I and case II are shown in table IV,V,VI,VII,VIII,IX,X and 
XI. The parameter values which lie outside the 3a range of global fit data are shown in red. 
The results clearly show that the predictions for the parameters do not change substantially 
if we go from type I dominating case to type II dominating one. The same observation 
follows from the figures 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10 where we have shown the variation 
of the neutrino oscillation parameters with respect to 7 G (0.001,1.0) for all the cases of 
interest. The variation of the oscillation parameters with respect to 7 are almost identical 
for both type I dominating and type II dominating cases with one overlapping the other in 
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FIG. 2: Variation of the predicted values of Am^ as a function of 7 in NH case 

most of the cases as seen from the figures. 

Apart from the observation that predictions for neutrino oscillation parameters remain 
almost same for both type I and type II dominating cases, we also observe that certain 

combinations of IH or NH, Majorana CP phase pattern (H h) or (+ + +), uilr of type 

CL or UQ do not give rise to correct 3a predictions of oscillation parameters at 7 ~ 1 where 
according to our selection criteria, the model is most likely to survive. Such combinations 
are 

• IH with Majorana CP phase pattern (+ + +) for CL type m^jj (Table V). 

• IH with Majorana CP phase pattern (H h) for UQ type thlr (Table VI). 

• IH with Majorana CP phase pattern (+ + +) for UQ type (Table VII). 

• NH with Majorana CP phase pattern (H h) for UQ type ttilr (Table X). 

• NH with Majorana CP phase pattern (+ + +) for UQ type m LR (Table XI). 

Thus, the UQ type texture of Dirac neutrino mass matrix is disfavored in our analysis 
for both types of hierarchies and Majorana CP phase patterns. For CL type Dirac neutrino 
mass matrix, normal hierarchy is compatible with both types of Majorana CP phase patterns 
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FIG. 3: Variation of the predicted values of sin 2 612 as a function of 7 in NH case 



whereas inverted hierarchy gives correct predictions only for the (H h) pattern of Majorana 

CP phase. 

Also, for these three favored models IH(H — h), NH(H — h) and NH(+++) with CL type 
rriLRi the predictions for Am 2 3 (rH), Amg 1 (NH) remain very close to the best fit central 
values whereas IH prefers Am| x to be at the upper end of the 3a range. Regarding the 
mixing angles, sin 2 #23, sin 2 # 13 stay very close to the best fit central value for all the models 
whereas sin 2 # 12 stays at the lower end of 3a range for NH(+++) model and upper end of 
3a range for IH(+-+), NH(+-+) models. Therefore, it is very likely that more precise future 
data from neutrino oscillation experiments will rule out some (if not all) of the scenarios 
discussed in this work. 



IV. CONCLUSION 

We have studied the survivability of quasi-degenerate neutrino mass scenario within the 
framework of generic left-right symmetric models by considering both types of hierarchies 

(normal and inverted), two types of extremal Majorana CP phases (H h) and (+ + +), 

and two types of Dirac neutrino mass matrices: charged lepton type and up quark type. 
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FIG. 4: Variation of the predicted values of sin 2 #23 as a function of 7 in NH case 
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In generic LRSM, neutrino mass can get contributions from both type I as well as type 
II seesaw terms. We check that for generic choices of left-right symmetry breaking scale 
and Dirac neutrino mass matrices, these two different contributions to neutrino mass can 
be comparable only if a dimensionless parameter appearing in the type II term (7 in our 
notation) is fine tuned to be very small. Without considering this special fine-tuned case, 
we consider two other possible cases: one in which type I seesaw dominates and the other 
in which type II seesaw dominates. We use the generic parametrization of neutrino mass 
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FIG. 10: Variation of the predicted values of sin 2 #13 as a function of 7 in IH case 



matrix(TBM plus non leading corrections giving rise to non-zero #13) obtained by several 
groups and find the numerical value of these parameters using the global fit neutrino oscilla- 
tion data as well as the cosmological upper bound on the sum of absolute neutrino masses. 
We then keep the dominating seesaw term fixed in the seesaw formula and vary the other 
term by varying the dimensionless parameter 7. We then calculate the predictions for neu- 
trino oscillation parameters for different values of 7 and check whether they agree with the 
3o" range of global fit data at 7 ~ 1 which corresponds to the case where the non-leading 
term in the seesaw formula has the minimum possible contribution. 

Apart from the observation that both type I dominating and type II dominating cases 
give almost identical predictions, we also observe that up quark type Dirac neutrino mass 
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matrix is disfavored within our framework. Predictions for neutrino parameters using up 

quark type with both IH, NH and (H h), (+ + +) deviate from the 3a range of 

global fit data. In the case of charged lepton type m^g, normal hierarchy survives with 
both types of Majorana CP phase patterns whereas inverted hierarchy survives only with 

the (H h) case. In view of above, the scenario of quasi-degenerate neutrinos can survive in 

nature within the framework of type I and type II seesaw mechanism and hence can not be 
ruled out yet apart from certain exceptions mentioned above. However, some of the model 
predictions lie at the extreme ends of the 3a allowed range and hence more precise data 
from neutrino oscillation experiments should be able to to rule out some (if not all) of the 
scenarios we have discussed in our work. 
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